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I. INTRODUCTION 



In recent years there has been extensive research on 
magnetic molecules as these are novel, realizable systems 
for exploring magnetic phenomena in low-dimensional 
magnetic materials J^'^'^'^ Among these diverse systems 
the pair of Keplerate structural type magnetic molecules 
abbreviated as {Mo72Fe3o} and {Mo72Cr3o}, each host- 
ing a highly symmetric array of 30 exchange-coupled 
magnetic ions ("spin centers"), serve as highly attrac- 
tive targets for the investigation of frustrated magnetic 
systems. In these molecules-'^ the magnetic ions Fe^'^^ 
(spin s = 5/2) and Cr^^^ (spin s = 3/2) occupy the 
30 symmetric sites of an icosidodecahedron, a closed 
spherical structure consisting of 20 corner-sharing trian- 
gles arranged around 12 pentagons (diamagnetic polyox- 
omolybdate fragments). This is a zero-dimensional ana- 
logue of the planar kagome lattice that is composed of 
corner-sharing triangles arranged around hexagons. A 
useful theoretical framework that has been employed^i^iiS 
in studying these magnetic molecules is based on an 
isotropic Heisenberg model, where each magnetic ion is 
coupled via intra-molecular isotropic antiferromagnetic 
exchange to its four nearest neighbor magnetic ions, and 
all of the 60 intra-molecular exchange interactions are 



of equal strength (henceforth, "single- J model" )jii Un- 
fortunately the quantum Heisenberg model of the two 
magnetic molecules is intractable using either analyti- 
cal or matrix diagonalization methods. Nevertheless, the 
nearest-neighbor exchange constant for each molecule has 
been established by comparing experimental data above 
30 K for the temperature-dependent zero-field suscepti- 
bility with data obtained by simulational methods us- 
ing the quantum^ and classical^ Monte Carlo methods. 
For temperatures below about 30 K, where the quantum 
Monte Carlo method proves to be ineffective for the two 
magnetic molecules due to frustration effects the clas- 
sical Heisenberg model is at present the only practical 
platform for establishing the dependence of the magne- 
tization M [H, T) on external magnetic field H and tem- 
perature T. A rigorous analytical result^ for the classical, 
nearest-neighbor, single-J Heisenberg model states that, 
in the zero temperature limit, M is linear in H until M 
saturates (saturation fields Hs — 17.7 T and 60.0 T for 
{Mo72Fe3o} and {Mo72Cr3o}, respectively). The practi- 
cal relevance of the classical Heisenberg model in describ- 
ing these magnetic molecules even at low temperatures 
was strikingly demonstrated in an earlier experiment^ on 
{Mo72Fe3o} at 0.4 K, showing an overall hnear depen- 
dence of M on and its saturation at approximately 
17.7 T. The ground state envisaged by the classical single- 
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J model is characterized by high-symmetry spin frustra- 
tion. In particular, for H — 0, in the ground state the 
spins are coplanar with an angular separation of 120° 
between the orientations of nearest-neighbor spins. On 
increasing the external field H the spin vectors gradu- 
ally tilt towards the field vector, until full alignment is 
achieved when H — H^, while their projections in the 
plane perpendicular to the field vector retain the 120° 
pattern for nearest-neighbor spins4 It would be reason- 
able to expect that M{H, T) is an analytic function of 
its variables and thus dM/dH w Ms/Hs for H < 
and ksT <C Jo, essentially independent of both H and 
T in these intervals, where Jo is the exchange constant 
of the single- J classical Heisenberg modelj^ Indeed our 
classical Monte Carlo calculations confirm this behavior. 

As shown here, new and crucial features of M{H,T) 
become accessible upon examining the differential sus- 
ceptibility, dM/dH, as a function of H and T. What 
emerge are significant conflicts, spelled out in Sees. Ill Al 
and III Bl between the results of our measurements and 
a theory based on the classical single-J model as con- 
cerns both the T and H dependence of dM/dH below 
5 K. However, and this is the central idea of this pa- 
per, agreement with experiment is achieved upon adopt- 
ing a refined classical Heisenberg model of {Mo72Fe3o} 
and {Mo72Cr3o}, where we drop the assumption of 
a single, common value for the nearest-neighbor ex- 
change constant, using instead a rectangular probability 
distribution's, for the 60 nearest-neighbor intra-molecular 
interactions. 

The layout of this paper is as follows. In Sec. Ill Al we 
present our experimental results for dM/dH as a func- 
tion of T in the low-temperature range for {Mo72Fe3o} 
and {Mo72Cr3o} for several low field values. In Sec. Ill Bl 
we present our experimental results for dM/dH as a 
function of H. In Sec. IIIII we show that the model 
system of a classical isosceles spin triangle exhibits the 
same qualitative features as the results of our experi- 
ments summarized in Sec. Ill Al In Sec. IIVI using the 
classical Heisenberg model of the icosidodecahedron and 
Monte Carlo simulational methods, we show that excel- 
lent agreement between a model calculation based on 
a two-parameter distribution of exchange constants and 
our experimental data of Sec. |TT] can be achieved. Fi- 
nally, in Sec. |V] we summarize our findings, discuss the 
broader implications of our results, and identify several 
open questions. 



II. EXPERIMENT 
A. dM/dH versus T 



paramagnetic impurities. Our data was obtained by 
using a self-resonating LC circuit driven by a tunnel 
diodej ^^i^^ Briefly, a tank circuit consisting of a small coil 
of inductance Lq and a capacitor, C, is kept at constant 
temperature of 5 ± 0.005 K. The sample is mounted on 
a sapphire holder, which is inserted into the coil without 
making contact. In the absence of the sample, the cir- 
cuit resonates at the resonant frequency 27r/o = 1/ ^ LC . 
When a sample with susceptibility x is inserted into the 
coil, the resonant frequency changes from /q to /(x) due 
to the change of the coil inductance, L ~ d^/dl, where 
$ is the total magnetic flux through the coil and / is the 
current in the coil. This current generates a magnetic 
field of about 20 mOe. If the magnetic perturbation due 
to the sample is small, the shift of the resonant frequency 
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Measurements of dM/dH versus T were performed 
at Ames Laboratory on polycrystalline samples of 
{Mo72Fe3o} and {Mo72Cr3o} (see experimental sections 
of Refs. 'a and ^ prepared by employing optimized syn- 
thesis methods and re-crystallization steps to minimize 



FIG. 1: Temperature dependence of the measured differen- 
tial susceptibility dM/dH for {Mo72Fe3o} (upper panel) and 
for {Mo72Cr3o} (lower panel). Values of H are listed in the 
legends. 
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Here N is the demagnetization factor, AL (x) = 
|L (x) — io| Lq is the change of the coil inductance, 
V is the sample volume and Vb is volume of the coil, 
and X is the dynamic magnetic susceptibility at our typ- 
ical resonant frequency of 10 MHz. This frequency is 
still much lower than characteristic frequencies and the 
response can be considered as static, i.e., x can be iden- 
tified with dM/dH . Measured frequency shifts for the 
samples described in this work are of the order of 1 to 
10 Hz, whereas the experimental resolution of the setup is 
about 0.01 Hz which corresponds to a smallest detectable 
magnetic moment of about 5 picoemu. For easy compari- 
son with other experiments, the measured frequency shift 
is proportional to the change in the total magnetic mo- 
ment of the sample at a given frequency, AM = CA/, 
where C is a calibration constant. 

In Fig. [T] we show our experimental data for dM/dH 
versus T for {Mo72Fe3o} and {Mo72Cr3o}, respectively. 
The calibration of the experimental data was achieved by 
matching the effective magnetic moment inferred from 
the measured frequency shift to low-field (3 - 50000 
Gs) ac (10 to 1000 Hz) and dc susceptibility measure- 
ments performed on a Quantum Design MPMS, which 
show that at high temperatures dynamic effects can be 
neglected. The assumption of static behavior is sup- 
ported by the fact that all curves collapse onto each 
other for T > 5 K. The most striking feature of these 
data is the very strong dependence of dM / dH on H 
and T. As remarked in Sec. [H this behavior is con- 
trary to that predicted by the single- J model, specifically 
M « Hxo{kBT/Jo) « Hxo{0), that is, dM/dH is inde- 
pendent of H and T in the regime H Hg, ksT <^ Jq. 



B. dM/dH versus H 

In Fig. [3] we show our experimental data for dM /dH 
versus H for both {Mo72Fe3o} (upper panel) and 
{Mo72Cr3o} (lower panel). These data sets were obtained 
by measuring the magnetization using pulsed magnetic 
fields (typical sweep rate 15000 T / s) with a standard in- 
ductive method at facilities of Okayama University, To- 
hoku University, and University of Tokyo. Also shown 
are the corresponding classical Monte Carlo simulational 
results for the single- J model (dashed curve) for the in- 
dicated temperatures. 

The inadequacy of the single-J model is striking in 
that the simulational data differ from the experimental 
data in four important ways. First, the experimental 
data exhibit a steep rise for decreasing low fields, and the 
lower the temperature the steeper the rise. This is con- 
sistent with our experimental findings in Sec. Ill Al for the 
T- dependence of dM/dH. These features are entirely 
absent for the single-J model; in particular, dM/dH is 



essentially independent of field in the low-field regime. 
Second, the local minimum in the experimental data is 
significantly broader than that predicted by the single-J 
model. Third, according to the single-J model a local 
minimum in dM/dH versus Jfemcrges at T = K at a 
field H = H,/3 (approx. 6 T and 20 T for {Mo72Fe3o} 
and {Cr72Fe3o}, respectively) and it is enhanced with 
increasing T.^-^ By contrast, our experimental data for 
fields in the vicinity of iJs/3 differ insignificantly with 
temperature. Fourth, for {Mo72Fe3o} the experimental 
data for 0.42 K shows a decrease with increasing field 
above 10 T, quite distinct from the pattern of the single-J 
model. 



III. CLASSICAL ISOSCELES SPIN TRIANGLES 

In this Section we give a qualitative explanation for 
our experimental findings in Sec. Ill Al We suggest that 
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FIG. 2: Magnetic field dependence of the measured differen- 
tial susceptibility dM/dH, shown in red, for {Mo72Fe3o} for T 
= 0.42 K and 60 mK (upper panel, inset) and for {Mo72Cr3o} 
for T = 1.3 K and 0.5 K (lower panel, inset). The dashed 
curves are the results of the single-J model for these temper- 
atures. 
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the strong sensitivity of the differential susceptibihty on 
H and T reflects a non-analytic dependence of the mag- 
netization on these variables, a characteristic already ex- 
hibited by independent classical isosceles spin triangles. 
The Hamiltonian of a spin triangle is given in footnote l3 
in terms of dimensionlcss quantities. In particular, the 
three pair- wise interactions are described by two different 
antiferromagnetic exchange constants (positive values of 
J and J') and we consider both cases, J' < J and J' > J. 

For the equilateral spin triangle {J' = J — Jq) at T = 
K the magnetic moment per triangle, M{H,0), is a 
continuous, linear function of H and, in particular, it 
vanishes for H ^ 0. Specifically, M{H, 0) = 3H/Hs for 
—Hs < H < Hs, where Hs = 3 Jo is the saturation field. 

Analytical calculation of M{H, 0) for the correspond- 
ing isosceles spin triangle is a non-trivial task for general 
values of J/ J' as it is necessary to carefully identify the 
configuration of least energy for arbitrary values of H . 
The final results are as follows: For the case J' > J: 
M{H > 0, 0) = 1 - J/J' + (2 + J/J')H/Hs for < JT < 
Hg, where Hg — 2 J' + J is the saturation field. For the 
case J' < J < 2 J': M{H > 0, 0) = J/ J' - 1 + H/J' for 
< iJ < 2 J' - J; M{H, 0) = 1 for 2 J' - J < Jf < J; 
and M{H,0) = H/J for J < iJ < 3J = iJ,. M{H,Q) 
is subject to the antisymmetry property M{—H, 0) = 
—M{H,0). The discontinuous behavior for Jf ^ is of 
particular importance in both cases. Although this dis- 
continuity occurs only at K, it is the crucial mechanism 
responsible for the great sensitivity of AI {H, T) in the low 
H - low T regime. 




FIG. 3: Temperature dependence of the differential suscep- 
tibility dM/dH for a classical isosceles spin triangle (black 
curves) with J/fcs = 1 K and J' jkn = 1.8 K for field val- 
ues = 0, 0.1, . . . , 1. Inset: MjMs versus H for the 
values fcsT/Jo = 0,0.001,0.01,0.02,0.05,0.1,0.2,0.5. Re- 
sults for the corresponding classical equilateral spin triangle 
{Jo/ks = 1.267) K are given by the red dashed curves. 

The form of M{T, H) for finite T can, in principle. 



be derived for this model system by analytical methods, 
however these calculations are substantially more intri- 
cate than for the classical equilateral spin triangle (see 
SeciniB of O. Ciftja et al. listed in Ref.Ull). For practi- 
cal purposes, the simplest procedure is to use the classi- 
cal Monte Carlo method for convenient numerical choices 
of J'/ J. In Fig. Owe display our results for the choices 
J/fcs = 1 K, J' Iks = 1-8 K. The most important feature 
to note is that while M is indeed a continuous function of 
iJ for any nonzero T, the quantity dM/dH, provided in 
the main portion of Fig. [31 exhibits strong temperature 
dependence for weak magnetic fields. Indeed, the curves 
fan with increasing H m. a manner that is strikingly sim- 
ilar to that shown in Fig. [T] Similar behavior occurs for 
cases where J' < J. By contrast, for the correspond- 
ing classical equilateral spin triangle dM / dH is virtually 
independent of T, as is seen in Fig. [S] 



IV. MULTIPLE NEAREST-NEIGHBOR 
COUPLINGS 

The T and H dependence of dM/dH for the classical 
isosceles spin triangle considered in the previous section 
is remarkably similar to that of our experimental data 
in Sec. ini However, to achieve a more realistic model, 
in the following we assume that the 60 nearest-neighbor 
exchange interactions between magnetic ions in a given 
molecule are characterized by a probability distribution 
with two adjustable width parameters. In the following 
Section we rationalize the use of a probability distribu- 
tion as a convenient way for summarizing the combined 
effects of multiple microscopic mechanisms that disrupt 
the use of an idealized, single- J model. 

We simulate each of {Mo72Fe3o} and {Mo72Cr3o} 
by considering an ensemble of up to 100 indepen- 
dent systems, for a total of 60 couplings per system. 
We assign values of the 6000 exchange constants us- 
ing a random number generator according to the fol- 
lowing rules: 1) The average value, Jo„, of the clas- 
sical exchange constant (in units of Boltzmann's con- 
stant) for the nth system is allowed to assume any value 
in the interval ((1 — r) Jq, (1-f t)Jo) with equal prob- 
ability, where Jq is chosen as 13.74 K for {Mo72Fe3o} 
and 32.63 K for {Mo72Cr3o} as determined by high- 
temperature susceptibility measurements^^; 2) For the 
nth system, the individual values of the classical ex- 
change constant are allowed to assume any value in the 
interval ((1 — p) Jon, (1 + p)Jon) with equal probability. 
For each molecule the two parameters r, p characteriz- 
ing these rectangular probability distributions were de- 
termined so as to provide an optimal fit with our exper- 
imental data for dM /dH versus H . 

In Fig. m we present our results for dM/dH versus 
H for {Mo72Fe3o} and {Mo72Cr3o}. Note the excellent 
agreement between the experimental data and the sim- 
ulational results obtained using our multiple-J model 
(solid curve). In the case of {Mo72Fe3o}, the opti- 
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mal choices of the parameters r, p were r = 0.15 and 
p = 0.40, whereas for {Mo72Cr3o} these were t — Q and 
p = 0.5. 

The motivation for using two distributions as described 
above is the foUowing. One can attribute the influence of 
the two distributions, characterized by r and p, respec- 
tively, to complementary effects which only in combina- 
tion lead to the observed properties of both molecules. 
The value of r controls the variation in the values of 
the mean exchange constant per molecule^ which leads 
to variations in the value of the saturation field Hg and 
hence in the value of the minimum in dM/dH versus 
H at Hs/S. By averaging over those variations one finds 
that dM/dH versus H starts to decrease at a much lower 
value of H than predicted by the single- J model and si- 
multaneously finds that the minimum at Hs/S is broad- 
ened as observed. However, using this distribution alone 
one cannot explain the observed strong H dependence 
of dM/dH in the low T - low H regime, because each 



molecule is still characterized by a single exchange con- 
stant. Introducing a second distribution, characterized 
by the parameter p, leads to a variation in the values of 
the exchange constant within a molecule with the effect 
that the corner-sharing spin triangles are of the scalene- 
type rather than equilateral-type. This gives rise to the 
non-analytic behavior of the magnetization aX H — Q for 
K and hence to the characteristic effects we have found 
in the low H - low T regime. In our simulations we have 
studied a very large range of choices of parameter pairs. 
The optimal choice for {Mo72Fe3o} can be narrowed to 
r = 0.15 ± 0.02 and p = 0.40 ± 0.02. For {MoyaCrao} we 
find p = 0.50±0.02, however r can be chosen in the range 
to 0.2 without any observable effect. This is due to the 
fact that for {Mo72Cr3o} magnetization measurements 
above the saturation field {Hg = 60 T) are not achiev- 
able at the present time. In the case of {Mo72Fe3o}, for 
which Hs = 17.7 T, the availability of magnetization data 
above the saturation field allows for a greatly reduced un- 
certainty in the value of r. 
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FIG. 4: Measured differential susceptibility dM/dH versus 
H, shown in red, for {Mo72Fe3o} for T = 0.42 K and 60 mK 
(inset) and for {Mo72Cr3o} for T = 1.3 K (inset: 0.5 K) 
and simulational results (solid black curve) using a multiple- J 
model for the optimal choice of the probability distribution 
parameters as given in the text. 
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FIG. 5: Simulational results for dM/dH versus T based on 
the multiple- J model using the optimal distribution parame- 
ters for {Mo72Fe3o} and for {Mo72Cr3o} as given in the text. 
Results for the corresponding single- J model calculations are 
shown as red dashed curves. 



Shown in Fig. [5] are our simulational results for 
dM/dH versus T for several different values of H using 
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the probability distribution with optimal parameters ap- 
propriate for {Mo72Fe3o} and {Mo72Cr3o}. These results 
(black curves) are strikingly similar to the experimental 
curves seen in Fig. [TJ the corresponding curves (shown 
in red) for the single- J model, for the same choice of the 
mean value Jo , are essentially indistinguishable from one 
another. This again strongly supports the existence of 
the multiple-J scenario. Note that, in the experiment 
as well as in the simulations, with increasing tempera- 
ture the curves for different field values rapidly converge 
and become indistinguishable from one another. Also, 
for increasing temperature, the results for the multiple-J 
model merge with those of the single-J model, as ex- 
pected, since the average exchange constant across the 
ensemble, Jq, is chosen to equal to the exchange con- 
stant of the single-J model. Finally, it remains to be 
seen whether the sharp rise in the curves of the upper 
panel of Fig. [S] below 200 mK is an experimental feature 
in {Mo72Fe3o} or merely an artifact of the multiple-J 
model based on a rectangular probability distribution. 



V. SUMMARY AND DISCUSSION 

In this article, we have presented our experimental 
data for the differential susceptibility of the pair of mag- 
netic molecules {Mo72Fe3o} and {Mo72Cr3o} as a func- 
tion of magnetic field and temperature. Below 5 K these 
data are strikingly different from what can be provided 
using a classical Heisenberg model with a single value of 
the nearest-neighbor exchange constant (single-J model). 
We have achieved excellent agreement with our experi- 
mental data upon adopting a classical Heisenberg model 
where the 60 nearest-neighbor interactions are not identi- 
cal; instead, the values of the exchange constants are de- 
scribed by a two-parameter probability distribution with 
a mean value as determined from experimental dM/dH 
data above 30 K using the single-J model. Above 5 K 
the single-J model provides a satisfactory description of 
each molecule. 

Since the icosidodecahedron structure consists of 
corner-sharing triangles, it is not surprising that the 
Heisenberg model of independent classical isosceles spin 
triangles provides a simple yet instructive model in that 
it exhibits the main qualitative features of our experi- 
mental data. We note here that a similar approach has 
been employed successfully for various two-dimensional 
spin systems on triangular lattices^S. For example, in the 
case of manganese tricyanomethanide the so-called 'row 
model' based on connected isosceles triangles provides the 
explanation for an unusual magnetic-field dependence of 
the spin ordering^^. In the context of independent clas- 
sical isosceles spin triangles one can figuratively describe 
the effect of multiple exchange constants as modifying the 
spin frustration from the standard 120° angular separa- 
tion between spin vectors of the equilateral spin triangle. 
The operational consequence is that the magnetization, 
for T = 0, of an isosceles spin triangle is a non- analytic 



function of magnetic field for TJ = 0, and this is mani- 
fested in dM/ dH being a highly sensitive function of its 
arguments for small H and T . 

The existence of a distribution of nearest-neighbor ex- 
change constants can be expected to be responsible for a 
significant lifting of degeneracies of magnetic energy lev- 
els. To be specific, the quantum rotational band model^-^, 
which is a solvable alternative to the nearest-neighbor 
single-J quantum Heisenberg model, predicts a discrete 
spectrum of energy levels, many of which have a very 
high degeneracy due to large multiplicity factors. Per- 
turbing this model Hamiltonian by using a distribution 
of J-values would remove a major fraction of these de- 
generacies. The lifting of level degeneracies could provide 
a reasonable explanation for three long-standing puz- 
zling issues concerning these magnetic molecules: The 
first issue is the very broad peak (maximum at 0.6 meV) 
that has been observed by inelastic neutron scattering on 
{Mo72Fe3o} at 65 mK. In order to qualitatively repro- 
duce the observed peak using the rotational band model, 
it was necessary in Ref. [2l| to perform the calculations 
upon assigning a large energy width (0.3 meV) for the 
individual energy levels. We suggest that the source of 
this large energy width might be the lifting of the major- 
ity of degeneracies associated with a single-J model. 

A second important consequence of the splitting of 
highly degenerate levels would be that the molecules 
could exhibit classical characteristics down to very low 
temperatures. This would provide a very reasonable ex- 
planation for the surprising fact that our simulational re- 
sults based on the classical Heisenberg Hamiltonian are 
so successful in describing {Mo72Cr3o}, despite the fact 
that the Cr^'^^ ions have a small spin (3/2). Stated differ- 
ently, with the lifting of degeneracies and the fanning out 
of energy levels the effective temperature for the crossover 
from classical to quantum behavior can be anticipated to 
be considerably lower than that expected a priori for the 
single-J model. 

Third, the failure of efforts to observe magnetiza- 
tion steps, in measurements of magnetization versus iJ, 
in the mK temperature range in both {Mo72Fe3o} and 
{Mo72Cr3o} could also be attributed to the removal of de- 
generacies of the magnetic energy levels. The occurrence 
of magnetization steps at low temperatures is associated 
with the field-induced crossing of successive energy levels 
of the lowest rotational band. However, a discrete level 
associated with total spin quantum number S has multi- 
plicity 25+1 [total degeneracy (25' -1-1)^]. If this degener- 
acy is lifted there will be a multitude of level crossings at 
slightly different field values and thus give rise to blurred 
effects down to lower temperatures than would otherwise 
be expected. 

Given the finite-spin values of the Fe^^^ and Cr'^^^ ions, 
is it possible to explain the present experimental find- 
ings based on a quantum Heisenberg model that adopts 
a common, single value of the exchange constant for all 
of the nearest-neighbor interactions? We strongly doubt 
that this is possible, for we have seen, albeit with a clas- 



7 



sical Heisenberg model, that it is the spread in values 
of the nearest-neighbor exchange constant that fuels the 
sensitive dependence of dM/dH on T, or equivalently 
the non-analytic behavior of magnetization on H in the 
low H - low T regime. 

Basing our simulations on a probability distribution 
for the nearest-neighbor exchange interaction has led to 
excellent agreement with the detailed features of our ex- 
perimental data including the sensitive dependence of 
dM/dH on T and H. One can attribute the failure 
of the single- J model to the combined effect of a large 
number of diverse perturbing mechanisms. The effects 
of impurities, variations in the exchange-coupling geome- 
try, weak magnetic exchange interactions of more-distant 
neighbors, Dzyaloshinsky-Moriya and dipole-dipole in- 
teractions in these magnetic molecules are some of the 
many effects that are excluded when one uses an ide- 
alized single-J model. On the other hand, it is at this 
stage an extremely difficult, essentially impossible task 
to realistically quantify the effects of the diverse mecha- 
nisms. A theoretical description based on a Heisenberg 
model where the nearest-neighbor exchange constant is 
chosen using a probability distribution provides a rel- 
atively simple, phenomenological platform for compro- 
mising between the need for microscopic realism versus 
practical limitations. Ultimately it is significant that a 
two-parameter probability description can actually pro- 
vide the level of agreement that we have found. Finally, 
we remark that other choices of probability distributions 
can be expected to perform equally well. 

As one example of the complications in assessing the 
plethora of perturbing mechanisms, we consider the vari- 
ation in the intramolecular distances between nearest- 
neighbor magnetic ions. A geometric analysis utiliz- 
ing existing low-temperature single crystal X-ray struc- 
ture data for {Mo72Fe3o} molecules shows that the sub- 
structure of the magnetic ions is close to an ideal It- 
symmetric geometry. There is a standard deviation of 
0.04525 A (0.70%) and a maximum deviation of 1.4% 
from the average Fe-Fe distance of 6.4493 A for aU 60 
Fe-Fe nearest-neighbor distances. Besides the distances 
between the spin centers, geometric variations within the 
polyoxomolybdate exchange ligand are observed. In par- 
ticular the 0-Mo-O bond angle variations should affect 
the total orbital overlap and thus the exchange energy 
due to the spatially anisotropic character of the (unoccu- 
pied) Mo(4d) orbitals. For {Mo72Fe3o}, an angular range 
103.6° — 106.6° is observed, which in part is caused by 
crystallographic disorder of Mo positions. To assess the 
influence of various geometric parameters involved in the 
superexchange pathways between two nearest-neighbor 
spin centers, both binding to a pentagonal diamagnetic 
[Moy^02i(H20)6]^" = {Moe} fragment, we performed 
systematic Density Functional Theory-Broken Symme- 
try calculations^ on a model system in which two s 
= 1/2 [VlV0(H20)2]2+ groups are coordinated to such 
a {Mog} fragment in a nearest-neighbor (1,2) configu- 
ration. The fragment is augmented by an additional 



Zn (H20)4 group binding in a (1,3,5) configuration for 
charge neutrality. The geometry of this model system 
was adjusted to match the actual configurations occur- 
ring in {Mo72Fe3o}. We find that J in such a model sys- 
tem can deviate by up to ±8% from the average value Jq. 
Given the similarity between V0^+, Cr^^'^, and Fe^^'^ in 
the Keplerate systems, namely that the magnetic orbitals 
cause a nearly isotropic spin density distribution, we ex- 
pect that the variations in the relative values of J span a 
very similar interval for {Mo72Fe3o} and {Mo72Cr3o}. As 
the intra- molecular variation in the values of the nearest- 
neighbor exchange constants implied by the optimal val- 
ues (given in Sec. lIV[) of the parameter p are significantly 
larger, we suggest that this is due to numerous other per- 
turbing mechanisms, some of which we listed above. 

We also note that other attempts to explain limited 
features of dM/dH, specifically the broadening of the 
minimum versus H for {Mo72Fe3o}, have been considered 
in the literature. One attempt assumed an elevated spin 
temperature during the pulsed field measurements, how- 
ever this could be ruled out since a subsequent steady- 
field measurement reproduced the results obtained by 
the pulsed- field technique. Second, in a simulational 
study based on classical Monte Carlo calculations, ef- 
fects of magnetic anisotropies, Dzyaloshinsky-Moriya, 
and dipole-dipole interactions have been consideredi^ 
However, our own comprehensive simulational studies of 
these same mechanisms have shown that they give rise to 
only very minor corrections on the width of the minimum 
in dM/dH versus H for any reasonable choices of model 
parameters. 

Finally, we suggest the additional possibility that in 
these molecules the variation of the exchange inter- 
action is spontaneously generated so as to lower the 
system's magnetoelastic energy. Such behavior has 
been observed experimentally and described theoreti- 
cally for a variety of antiferromagnetic oxide pyrochlore 
compounds^^i^li^. The pyrochlore lattice consists of 
corner-sharing tetrahedra and exhibits geometric frus- 
tration. In this regard one can understand this struc- 
ture as the three-dimensional 'cousin' of the corner- 
sharing triangle type structures realized in {Mo72Fe3o} 
and {Mo72Cr3o}. 

In any event, it is highly satisfying that the frustrated 
magnetic molecules {Mo72Fe3o} and {Mo72Cr3o}, osten- 
sibly zero-dimensional systems, are a source of novel and 
intriguing magnetic behavior. 
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